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ORIGINAL ARTICLE 

Genome-wide association study identifies a potent locus 
associated with human opioid sensitivity 

D Nishizawa\ K Fukuda^ S Kasai\ J Hasegawa\ Y Aoki^'^ A Nishi\ N Saita^ Y Koukita^ M Nagashima^ R Katoh^ Y Satoh^ 

M Tagami^ S Higuchi^ H Ujike^'^^ N Ozaki^'^^ T Inada^'^^ N Iwata^'^^ I Sora^'^°'^^ M lyo^^'^^ N Kondo^^'^^ M-J Won^^'^^ N Naruse^'^'^^ 

K Uehara-Aoyama^^'^^ M Itokawa^^ M Koga^^ T Arinami^^ Y Kaneko^ M Hayashida^^ and K Ikeda^'^^ 

Opioids, such as morphine and fentanyl, are widely used as effective analgesics for the treatment of acute and chronic pain. In 
addition, the opioid system has a key role in the rewarding effects of morphine, ethanol, cocaine and various other drugs. Although 
opioid sensitivity is well known to vary widely among individual subjects, several candidate genetic polymorphisms reported so far 
are not sufficient for fully understanding the wide range of interindividual differences in human opioid sensitivity. By conducting a 
multistage genome-wide association study (GWAS) in healthy subjects, we found that genetic polymorphisms within a linkage 
disequilibrium block that spans 2q33.3-2q34 were strongly associated with the requirements for postoperative opioid analgesics 
after painful cosmetic surgery. The C allele of the best candidate single-nucleotide polymorphism (SNP), rs2952768, was associated 
with more analgesic requirements, and consistent results were obtained in patients who underwent abdominal surgery. In addition, 
carriers of the C allele in this SNP exhibited less vulnerability to severe drug dependence in patients with methamphetamine 
dependence, alcohol dependence, and eating disorders and a lower 'Reward Dependence' score on a personality questionnaire in 
healthy subjects. Furthermore, the C/C genotype of this SNP was significantly associated with the elevated expression of a 
neighboring gene, CREBl. These results show that SNPs in this locus are the most potent genetic factors associated with human 
opioid sensitivity known to date, affecting both the efficacy of opioid analgesics and liability to severe substance dependence. Our 
findings provide valuable information for the personalized treatment of pain and drug dependence. 
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INTRODUCTION 

The opioid system has important roles in both antinociception 
and reward. ^'^ Therefore, opioids, such as morphine and fentanyl, 
are widely used not only as effective analgesics for the treatment 
of acute and chronic pain but also as abused drugs. The opioid 
system is also involved in the rewarding effects of morphine,^ 
ethanol,"^ cocaine,^ and various other drugs^"^ or behaviors. 
However, opioid sensitivity is well known to vary widely among 
individual subjects,^ resulting in differences in the effectiveness of 
opioid analgesics and vulnerability to dependence on opioids and 
other drugs or behaviors. Individual differences may be 
attributable to both genetic and environmental factors,^ ° 
although the relative influence of each of these factors can be 
diverse. To date, several candidate genetic polymorphisms have 
been reported to be associated with opioid sensitivity in human 
studies.^ However, such polymorphisms have not sufficiently 
explained the wide range of interindividual variance observed in 



the sensitivity to opioid analgesics. A genome-wide approach has 
not yet been adopted to explore the best candidates, although 
this approach has been applied to other pharmacogenomics- 
related traits. Several genetic polymorphisms have been found to 
be associated with the sensitivity to pharmacotherapies.^^"^'' In 
this study, we sought to comprehensively identify genetic 
polymorphisms in the human genome that could greatly 
contribute to individual differences in opioid sensitivity by 
conducting a genome-wide association study (GWAS) of healthy 
subjects and further analyses. 



MATERIALS AND METHODS 

Subjects 

Enrolled in this multistage GWAS were 355 healthy patients who were 
scheduled to undergo cosmetic orthognathic surgery for mandibular 
prognathism at Tokyo Dental College Suidoubashi Hospital. The surgical 
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protocol and subsequent postoperative pain management were funda- 
mentally the same as a previous study^^ and are detailed in the 
Supplementary Information. 

The subjects recruited in the additional analysis to confirm the 
association between the rs2952768 single-nucleotide polymorphism 
(SNP) and postoperative opioid analgesic requirements were 112 patients 
who underwent major open abdominal surgery at several related hospitals. 
The surgical protocol and subsequent postoperative pain management 
were fundamentally the same as a previous study"'^^ and are detailed in 
the Supplementary Information. Enrolled in the study to investigate the 
contribution of the rs2952768 SNP to the symptoms of drug dependence 
or related personality traits were 203 patients with methamphetamine 
(METH) dependence with clinical data that included their multisubstance 
abuse status, 438 patients with alcohol dependence with clinical data that 
included the number of drugs used, 228 patients with eating disorders 
with clinical data including the presence or absence of other psychiatric 
disorders such as substance dependence and 500 healthy volunteer 
subjects with personality profile data from the temperament and character 
inventory (TCI).^^"^^ To examine the mRNA expression levels of the 
METTL21A {FAI\/m9A) and CREB1 genes, 100 post-mortem human brain 
specimens, from which DNA and RNA were extracted for experimental use, 
were obtained from the Stanley Medical Research Institute (SMRI; 
Bethesda, MD, USA) as samples independent of those in the association 
study with opioid sensitivity (SMRI samples). 

All of the individuals included in the study originated from Japan, 
with the exception of those from whom the SMRI samples were 
obtained, whose racial background was mostly European American (see 
Supplementary Information). 

The study protocol was approved by the Institutional Review Boards at 
the related hospitals, Tokyo Institute of Psychiatry (currently Tokyo 
Metropolitan Institute of Medical Science) and the ethics committee of 
each participating institute of the Japanese Genetics Initiative for Drug 
Abuse.^^'^^ All of the subjects provided informed, written consent for the 
genetic studies. The detailed demographic and clinical data of the subjects 
are provided in Supplementary Tables 1, 5-8 and 10. 

Genotyplng 

After total genomic DNA was extracted from whole-blood samples using 
standard procedures, whole-genome genotyping was performed using the 
Infinium assay II with an iScan system (lllumina, San Diego, CA, USA) 
according to the manufacturer's instructions. The data for the whole- 
genome genotyped samples were analyzed using BeadStudio or 
GenomeStudio with the Genotyping module v3.3.7 (lllumina) to evaluate 
the quality of the results. In the data-cleaning process, the samples with a 
genotype call rate of <0.95 were excluded from further analyses. As 
a result, one sample was excluded from further analyses. Markers with a 
genotype call frequency of <0.95 or 'Cluster sep' (that is, an index of 
genotype cluster separation) of <0.1 were excluded from the subsequent 
association study. A total of 295 036 SNP markers survived the filtration 
process and were used for the GWAS (Supplementary Figure SI). 

For additional genotyping of the rs2952768 and rs2254137 SNPs, the 
TaqMan allelic discrimination assay (Life Technologies, Carlsbad, CA, USA) 
was mostly conducted after total genomic DNA was extracted from 
whole-blood or oral mucosa samples using standard procedures. For 
samples that were not appropriately genotyped by this assay, direct 
sequencing was alternatively adopted to genotype the rs2952768 SNP. A 
total of 1 12, 203, 438, 228, 500 and 105 DNA samples from patients who 
underwent major abdominal surgery, patients with METH dependence/ 
psychosis, patients with alcohol dependence, patients with eating 
disorders, healthy volunteer subjects and SMRI, respectively, were used 
to genotype the rs2952768 SNP. In addition, a total of 105 DNA samples 
from the post-mortem specimens for the expression analysis were used to 
genotype the rs2254137 SNP, although genotyping this SNP for other 
samples was not conducted because of the strength of the linkage 
disequilibrium (LD) with the rs2952768 SNP. The genotype distribution of 
the rs2952768 SNP in patients with METH dependence/psychosis, patients 
with alcohol dependence and patients with eating disorders is provided in 
Supplementary Table 9. 

Quantitative PGR procedure 

The SMRI RNA samples were treated with DNase I using the RNase-Free 
DNase Set (Qiagen, Hilden, Germany), and clean-up was then performed 
using the RNeasy MinElute Cleanup Kit (Qiagen). First-strand complemen- 
tary DNA for use in the real-time quantitative PCR was synthesized with the 



SuperScriptlll First-Strand synthesis system for quantitative reverse 
transcriptase-PCR (Life Technologies) with lOOng purified total RNA 
according to the manufacturer's protocol and diluted properly with 
diethylpyrocarbonate-treated H2O before the experiments. 

To perform real-time quantitative PCR with a LightCycler 480 (Roche 
Diagnostics, Basel, Switzerland), TaqMan Gene Expression Assays (Life 
Technologies) were used as a probe/primer set specified for the FAMIIQA 
{/\/IETTL21A) gene and CREBl gene and a probe/primer set for the ACTB 
gene, a house-keeping gene that encodes j6-actin. The expression level of 
the FAIVIIWA {METTL21A) gene or CREBl gene was normalized to the 
expression level of the ACTB gene for each sample, and relative mRNA 
expression levels were compared between the genotype subgroups for 
each gene. The experiments were performed in triplicate (separate 
experiments) for each sample, and average values were calculated for 
normalized expression levels. 



Statistical analysis 

A three-stage GWAS was conducted for the patients who underwent 
painful cosmetic surgery to investigate the association between opioid 
sensitivity and the 295 036 SNPs that passed the quality control criteria in a 
total of 353 subjects (118, 117 and 118 subjects for the first-, second- and 
final-stage analyses, respectively). As an index of opioid sensitivity, 
postoperative patient-controlled analgesia fentanyl use during the first 
24-h postoperative period was used because analgesic requirements likely 
reflect the efficacy of fentanyl in each individual. To explore the association 
between the SNPs and phenotype, linear regression analyses were 
conducted in each stage of the analysis, in which postoperative fentanyl 
use (|ig kg"\- log transformed) and the genotype data of each SNP were 
incorporated as dependent and independent variables, respectively. 
Additive, dominant and recessive genetic models were used for the 
analyses because of the previously insufficient knowledge about the 
genetic factors associated with opioid sensitivity. The GWAS procedure is 
summarized in Supplementary Figure SI. In the first-stage analysis of 118 
subjects, the SNPs that showed statistical P-values of <0.05 were selected 
as the candidate SNPs for the second-stage analysis among the 295 036 
SNPs. For these SNPs, the second-stage analysis was conducted; again, the 
SNPs that showed P<0.05 were considered potent candidates and 
selected for further final-stage analysis. In the final stage of the three 
stages, the Q-values of the false discovery rate were calculated to correct 
for multiple testing in addition to the P-values based on previous 
reports.^"^'^^ The SNPs that showed Q<0.05 in the analysis were considered 
genome-wide significant. 

Additional analyses were conducted using the samples of the patients 
who underwent major abdominal surgery, patients with METH depen- 
dence/psychosis, patients with alcohol dependence, patients with eating 
disorders, healthy volunteers with personality profile data and post- 
mortem specimens for the expression analysis. For all of the genotype 
data used in these analyses, the distributions were checked using the 
test, and the absence of significant deviation from the theoretical 
distribution expected from Hardy-Weinberg equilibrium was confirmed 
(Supplementary Table 13). 

In the analysis of the patients who underwent major abdominal surgery, 
the calculated total dose of rescue analgesics administered during the first 
24-h postoperative period was used as an index of opioid sensitivity. To 
explore the association between the SNPs and phenotype. Student's f-test 
or analysis of variance was performed, in which Bonferroni correction for 
multiple comparisons was used for the post hoc tests. For these analyses, 
postoperative analgesic use (jig kg"\- log transformed) and the genotype 
data of each SNP were incorporated as dependent and independent 
variables, respectively, Tests were performed to investigate the 
contribution of the SNPs to the vulnerability to the presence of serious 
symptoms of substance dependence. For the analyses in which the 
number of subjects in a cell in 2 x 2 contingency tables was <5, Fisher's 
exact tests were conducted instead of tests. In the analysis of healthy 
volunteers with personality profile data, raw TCI scores were properly 
processed (see Supplementary Information). To explore the association 
between the SNPs and phenotype, linear regression analyses were 
conducted, in which the endpoint TCI score (log transformed) and 
genotype data of each SNP were incorporated as dependent and 
independent variables, respectively. The correction of multiple tests for 
the analyses of the seven phenotypes was not performed in the present 
exploratory study. In the analysis of the post-mortem specimens for the 
expression analysis, the calculated expression level of the FAMllQA 
{METTL21A) gene or CREBl gene normalized to the ACTB gene for each 
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sample was used. To explore the association between the SNPs and 
phenotype, Student's f-test or analysis of variance was performed, in which 
Bonferroni correction for multiple comparisons was used as the post hoc 
test. For these analyses, the relative expression level and genotype data of 
each SNP were incorporated as dependent and independent variables, 
respectively. 

For all of the statistical analyses described above, SPSS 18.0J for 
Windows (International Business Machines Corporation, Armonk, NY, USA), 
gPLINK V. 2.050, PUNK v. 1.07 (http://pngu.mgh.harvard.edu/purcell/plink/; 
accessed 1 March 2012),^^ and Haploview v. 4.1^^ were used. The criterion 
for significance was set at P<0.05, with the exception of the GWAS. 
Statistical corrections for multiple tests, such as Bonferroni adjustments on 
the multiple parameters analyzed, were not performed in the present 
exploratory study after the GWAS because it would be too conservative for 
genetic association studies,^^ meaning that the likelihood of type II errors 
is increased by Bonferroni adjustments, and truly important differences 
would be deemed nonsignificant.^^ 

Additional in silico analysis 

Additional in silico analyses and the Internet links for the websites of the 
databases referenced are provided in the Supplementary Information. 



RESULTS 

Identification of a potent locus associated with human opioid 
sensitivity by GWAS in subjects who underwent painful cosmetic 
surgery 

We first explored the association between genetic variations and 
opioid sensitivity in a total of 353 healthy subjects who were 
scheduled to undergo cosmetic orthognathic surgery (mandibular 
sagittal split ramus osteotomy) for mandibular prognathism that 
involved the administration of opioid analgesics (Supplementary 
Table SI), in which the surgical procedure was uniform and thus 
the invasiveness and resultant pain would be regarded as 
homogeneous among the subjects. A GWAS was conducted as a 
consecutive three-stage analysis to identify potent SNPs asso- 
ciated with the requirements for an opioid analgesic, fentanyl 
(|igkg"^), during the 24-h postoperative period (Supplementary 
Figure S1). Consequently, 9, 12 and 10 SNPs were selected as the 
top candidates for additive, dominant and recessive models for 
each minor allele, respectively, after the final stage 
(Supplementary Tables S2-4). Among these, several SNPs mapped 
to 2q33.3-2q34 showed significant associations after the final 
stage with 24-h postoperative fentanyl requirements in the 
additive and recessive models (additive model: combined 
j5 = 0.293, nominal P = 8.044 x 10"^; recessive model: combined 
j5 = 0.553, nominal P = 9.382 x 10~^; Supplementary Tables S2-4). 
The observed P-values of these SNPs, calculated as -loglO 
(P-value), obviously deviated from the expected values from the 
null hypothesis of uniform distribution in the quantile-quantile 
plot for the entire sample (Supplementary Figure S2). The genes 
located in this region were found to include METTL21A (FAI\/m9A) 
and CREB1, encoding methyltransferase like 21 A and cyclic 
adenosine 3',5'-monophosphate (cAMP) responsive element- 
binding protein 1 (CREB1), respectively. 

For further fine mapping of this region, we used the remainder 
of the genotyped data and imputed the genotype data from the 
SNPs that spanned 300 kbp, including both the METTL21A and 
CREBl genes (Figure la), and analyzed the association. We then 
identified several other SNPs around the best candidate SNP, 
rs2952768, with nominal P-values from 1x10"^ to 1x10""^ 
(Figure lb). After the LD analysis, all of these SNPs were found to 
be included in an LD block that spanned approximately 134 kbp 
and were in strong LD {D' = ^y>O.S) with rs2952768 (Figure 1c). 
Total fentanyl use during the 24-h postoperative period was 
1 .080 ± 0.053, 1 .1 46 ± 0.049 and 1 .666 ± 0.083 |ig kg"^ (log trans- 
formed; mean ± s.e.m.) in subjects with the T/T, T/C and C/C 
genotypes, respectively (Figure 2a). 



Association of rs2952768 SNP with sensitivity to opioid analgesics 
in patients who underwent major open abdominal surgery 

To examine whether the SNPs identified in our GWAS generally 
affect individual differences in opioid sensitivity, we attempted to 
confirm the association between the rs2952768 SNP and post- 
operative opioid requirements in another cohort that underwent a 
different surgical procedure. The subjects recruited in this study 
were 112 patients who underwent major open abdominal surgery 
under combined general and epidural anesthesia (Supplementary 
Table S5),^^ mostly gastrectomy for gastric cancer and colectomy 
for colorectal cancer, which involves different modes of invasion 
from the orthognathic surgery and might cause different pain 
modalities. Appropriate doses of analgesics, mainly opioids such 
as morphine, buprenorphine, pentazocine and pethidine, were 
administered as rescue analgesics at the discretion of the 
surgeons whenever the patients complained of significant 
postoperative pain during the postoperative period, and the 
total dose administered was estimated for the association analysis. 
As a result, a significant difference in postoperative analgesic 
requirements was found between the subjects with the combined 
T/T and T/C genotype and subjects with the C/C genotype in the 
rs2952768 SNP. Interestingly, the subjects with the C/C genotype 
required significantly more analgesics than the subjects with the 
combined T/T and T/C genotype in the rs2952768 SNP 
(fno= -2.340, P = 0.021), a pattern similar to the one observed 
in the subjects who underwent cosmetic orthognathic surgery. 
Total analgesic use, equipotent with systemic fentanyl, during the 
24-h postoperative period was 0.359 ±0.073, 0.397 ±0.068 and 
0.741 ± 0.1 69 |ig kg"^ (log transformed; mean ± s.e.m.) in the 
subjects with the T/T, T/C and C/C genotypes, respectively 
(Figure 2b). 



Association of rs2952768 SNP with severity of drug dependence in 
patients in several cohorts 

The results suggested that the subjects with the C/C genotype in 
the rs2952768 SNP required more analgesics than the subjects 
with the other genotypes, attributable to the decreased effective- 
ness of opioid analgesics in both cohorts. Given the fact that the 
opioid system is involved in both rewarding and analgesic effects, 
one could assume that decreased opioid sensitivity may reflect 
the decreased rewarding effects of various drugs or behaviors and 
less liability to serious dependence. To test this hypothesis, we 
investigated the contribution of the rs2952768 SNP to the 
vulnerability to substance dependence in additional subjects with 
METH dependence, alcohol dependence and eating disorders 
(Supplementary Tables S6-8). In the initial case-control analyses, 
no associations in the genotypic and allelic distribution of this SNP 
were found between the sulDjects with psychiatric disorders and 
corresponding control subjects (Supplementary Table S9). How- 
ever, a significant difference in the genotypic distribution was 
found between the absent and present subgroups of polydrug use 
among the patients with METH dependence (x^ = 3.979, 
p = 0.046). Indeed, fewer polydrug abusers carried the C/C 
genotype compared with monodrug users (Table 1). A similar 
result was found in the patients with alcohol dependence. A 
significant difference in the genotypic distribution was found 
between the absent and present subgroups of drug use among 
patients with alcohol dependence (x^ = 3.860, P = 0.049, and 
= 3.039, P = 0.097, in the additive and recessive models, 
respectively), and a lower proportion of drug abusers carried the 
C allele or C/C genotype compared with the alcoholics without 
drug abuse (Table 1). Furthermore, a significant difference in allelic 
distribution was found between the absent and present 
subgroups of comorbid dependence among the patients with 
eating disorders (x^ = 3.985, P = 0.046, and / = 4.488, P = 0.034, 
in the additive and dominant models, respectively), and a lower 
proportion of patients with drug dependence carried the C allele 
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Figure 1. Candidate locus possibly associated with hunnan opioid sensitivity, (a) Illustration of the genes in the genonnic region fronn position 
208000000 to 208 300 000 on chronnosonne 2 in the HapMap database (http://hapnnap.ncbi.nlnn.nih.gov/index.htnnl.ja; accessed 1 March 
201 2). (b) Fine nnapping of the candidate region after the innputation-based association analysis. The circle, square and triangle plots represent 
the results fronn the additive, donninant and recessive nnodels, respectively. The area between the dotted vertical lines represents the genonnic 
position fronn 208 070 000 to 208 240 000 on chronnosonne 2. (c) State of linkage disequilibriunn (LD) between the SNPs in the genonnic position 
fronn 208 070 000 to 208 240 000 on chronnosonne 2, based on the genotype data of the subjects who are derived fronn the Japanese 
population and underwent cosnnetic orthognathic surgery. Nunnbers in squares in which two SNPs face represent the percentage of the 
values calculated fronn the genotype data of the SNPs. Squares with asterisks represent = 1. Only the values >0.70 are highlighted. 



connpared with the patients without drug dependence (Table 1). 
Although nonsignificant, likely because of the snnall sannple size, a 
marginal difference was observed in allelic distribution between 
the absent (x^ = 3.780, P = 0.052) and present subgroups with a 
comorbid alcoholic state among the patients with eating 
disorders, and a lower proportion of patients with an alcoholic 



state carried the C allele compared with the patients without an 
alcoholic state (Table 1). Altogether, these results showed that 
carriers of the C allele among the patients with psychiatric 
disorders, especially the C/C genotype, tended not to abuse 
polydrugs and not have comorbid alcohol or drug dependence. 
Although the association was nonsignificant in the recessive 
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Abdominal surgery: 




rs2952768 genotype 



rs2952768 genotype 



Figure 2. Association analysis between opioid analgesic requirennents and the rs2952768 SNR (a) Total dose of analgesics adnninistered per 
body weight {\xgkg~^; log transfornned) during the 24-h postoperative period after cosnnetic orthognathic surgery (nnandibular sagittal split 
rannus osteotonny). (b) Total dose of analgesics adnninistered per body weight {[ig kg~\- log transfornned) during the 24-h postoperative period 
after nnajor open abdonninal surgery. *Q<0.05, greater dose of analgesic adnninistered in the C/C genotype connpared with the T/C and T/T 
genotypes with genonne-wide significance; ^P<0.05, greater dose of analgesic adnninistered in the C/C genotype connpared with the T/C and 
T/T genotypes with nonninal significance. The data are expressed as nnean ± s.e.nn. 



Table 1. Distribution of genotypes of the rs2952768 polymorphism and odds ratio between the patient subgroups based on clinical status 


n 


Genotype 




Allele frequency 












OR (95% CI) 




T/T 


T/C 


C/C 


T 


C 






P 




C allele 


C/C + T/C genotype 


C/C genotype 












Genotype 


Allele 


Dominant 


Recessive 


(vs T) 


(vs T/T) 


(vs T/T + T/C) 


Polydrug use 
In METH dependence/psychosis patients 
Absence 53 19 22 12 
Presence 141 58 67 16 


0.556 
0.644 


0.444 
0.356 


3.979 
0.137 


2.262 
0.133 


0.450 
0.503 


3.979 
0.046* 


0.71 
(0.45-1.11) 


0.80 
(0.42 - 1 .54) 


0.44 
(0.19-1.00) 


Drug use 
In alcohol dependence patients 
Absence 391 172 166 
Presence 45 25 18 


53 
2 


0.652 
0.756 


0.348 
0.244 


3.908 
0.142 


3.860 
0.049* 


2.180 
0.140 


3.039 
0.097§^ 


0.61 
(0.37-1.00) 


0.63 
(0.34-1.17) 


0.30 
(0.07-1.26) 


Drug dependence 
In eating disorder patients 
Absence 200 85 
Presence 21 14 


93 
6 


22 
1 


0.658 
0.810 


0.343 
0.190 


4.552 
0.103 


3.985 
0.046* 


4.488 
0.034* 


0.793 
0.705§ 


0.45 
(0.20-1.00) 


0.37 
(0.14-0.96) 


0.40 
(0.05-3.16) 


Alcoholic state 

In eating disorder patients 
Absence 151 61 
Presence 70 38 


72 
27 


18 
5 


0.642 
0.736 


0.358 
0.264 


3.995 
0.136 


3.780 
0.052^ 


3.731 
0.053^ 


1.171 
0.279 


0.65 
(0.41-1.01) 


0.57 
(0.32-1.01) 


0.57 
(0.20-1.60) 


Abbreviations: METH, methamphetamine; n, the number of samples; OR, odds ratio; 95% CI, 95% confidence interval. * 
Fisher's exact test was presented instead of that from test. 


P<0.05, ^0.05 ^P< 0.1, 


^P-value from 



model for the C allele in all of the analyses, possibly because of a 
lack of statistical power caused by the limited sample size, the 
present results suggest that carriers of the C allele in this SNP have 
less inclination to abuse drugs, consistent with our assumption 
that various drugs of abuse have decreased rewarding effects in 
subjects with decreased opioid sensitivity, making these subjects 
less liable to the expression of symptoms of serious dependence. 



Association of rs2952768 SNP with 'reward dependence' score on 
a personality questionnaire in healthy subjects 
Another interest is whether this SNP also affects personality traits 
related to the reward system in healthy people. To address this 



issue, we investigated the association between this SNP and data 
from the TCI, a personality profiling questionnaire, in healthy 
volunteers (Supplementary Table S10). Intriguingly, among the 
seven dimensions of the TCI, a significant association was found 
only for reward dependence (RD) (novelty seeking: - 0.009055, 
P = 0.2995; harm avoidance: ^5 = 0.00431 7, P = 0.6819; RD: 
P= -0.0175, P = 0.03265; persistence: p= -0.007554, P= 0.6295; 
self-directed ness: p= -0.01274, P = 0.1 709; cooperativeness: 
j5= -0.004122, P=0.5302; self-transcendence: j5 = 2.68E-06, 
P = 0.9998; Supplementary Figure S3). The RD value decreased as 
the copy number of the carried C allele increased among the 
subjects. The RD value on the TCI was previously shown to be 
positively correlated with activity of the caudate head,^° which has 
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Figure 3. Relative nnRNA expression level of the candidate genes between each genotype subgroup of the rs2952768 SNP in post-nnortenn 
brains, (a) Results for the METTL21A gene, (b) Results for the CREB1 gene. NS, no significant association between relative nnRNA expression and 
genotype subgroup (P^O.05); ^P<0.05, greater level of nnRNA expression in the C/C genotype connpared with the T/C and T/T genotypes with 
nonninal significance. The data are expressed as nnean ± s.e.nn. 



been shown to be associated with reward processing.^^ Although a 
future confirnnatory study is needed, the present data suggest the 
possibility that this SNP is one of the predisposing factors that 
partially contribute to the developnnental differentiation of 
personality traits related to RD. 

Association of rs2952768 SNP with mRNA expression level of the 
CREBl gene 

Considering the fact that the rs2952768 and all of the neighboring 
SNPs that showed sinnilarly strong associations with opioid 
sensitivity were included in the LD block that comprised the 
METTL21A {FAI\/m9A) and CREBl genes, a subsequent issue is the 
impact of these SNPs on the function of these genes. Our database 
search estimated that several of the SNPs could putatively cause 
functional alterations (Supplementary Table S11). To pursue this 
issue, we examined the mRNA expression levels of these genes 
using real-time quantitative PCR with RNA samples extracted from 
post-mortem subject specimens and compared the mRNA expres- 
sion levels between the genotype subgroups for the rs2952768 
SNP, which were determined by genotyping the DNA samples 
extracted from the corresponding subjects. Although no significant 
association was found in the relative mRNA expression level of the 
METTL21A {FAMIWA) gene between the genotype subgroups 
(^^2,97 = 0-372, P = 0.690), a significant association was found in the 
relative mRNA expression level of the CREBl gene between the 
combined T/T and T/C genotype subgroup and C/C subgroup 
(Figures 3a and b). These results showed that the C/C genotype of 
this SNP was significantly associated with elevated CREBl mRNA 
expression (tg8 = -2.561, P = 0.012). 



DISCUSSION 

By conducting a multistage GWAS for the first time in healthy 
subjects who were treated with opioid analgesics, we identified a 
potent locus for opioid sensitivity that encompasses an LD block 
that includes the most significant SNP, rs2952768. The C allele of 
this SNP, especially homozygotes of the C allele, was associated 
with more analgesic requirements, suggesting the possibility to 
classify the patients into groups of high responders and low 
responders or non-responders to the given opioid, which would 
presumably correspond to the T/T or T/C and C/C genotype 
groups, respectively. Surprisingly, this SNP was also found to be 
significantly associated with postoperative analgesic requirements 
in subjects who also received opioids but underwent different 
surgical procedures. Moreover, the C allele with likely less opioid 
sensitivity in this SNP was associated with a lower risk for serious 



symptoms of substance dependence and a lower RD score on the 
TCI in healthy subjects. These results show that SNPs in this locus 
are the most potent genetic factors associated with human opioid 
sensitivity known to date, affecting both the efficacy of opioid 
analgesics and liability to severe substance dependence. 

Opioids exert their effects by binding to opioid receptors (that 
is, G-protein-coupled receptors) and triggering signaling transmis- 
sion to several downstream effectors, including inhibition of 
adenylyl cyclase, activation of G-protein-activated inwardly 
rectifying potassium channels, and inhibition of voltage-gated 
Ca^^ channels.^^"^"^ Inhibition of adenylyl cyclase inhibits the 
production of cAMP, thus decreasing the active form of protein 
kinase A, phosphorylating CREB, and decreasing gene expression 
in the nucleus related to the action of analgesia and reward.^^'^^ 
Indeed, the administration of cAMP intracerebrally or 
intravenously antagonized morphine analgesia in nontolerant 
and tolerant mice.^'' Moreover, all of the major behavioral effects 
of morphine, including analgesia, tolerance, reward, and physical 
dependence and withdrawal symptoms, were attenuated in mice 
that lacked adenylyl cyclase 5.^^ Meanwhile, increased CREB 
function in the nucleus accumbens decreased the sensitivity to 
the rewarding effects of morphine and cocaine in animals, 
whereas decreased CREB function had the opposite effect.^^'^^'^^ 
Although no change in opioid-binding sites or morphine-induced 
analgesia was observed in CREBa5-deficient mice,'^° alterations in 
dose-dependent morphine-induced reward were also reported in 
CREBa6-deficient mice."^^ Altogether, these reports indicate the 
involvement of CREB and the cAMP pathway in the analgesic and 
rewarding effects of opioids. Higher mRNA expression levels of the 
CREBl gene in subjects with the C/C genotype in the rs2952768 
SNP identified in our GWAS may indicate elevated CREB function 
and decreased sensitivity to the rewarding effects of opioids, 
resulting in greater postoperative opioid analgesic requirements 
and less vulnerability to dependence on other drugs. 

The best candidate SNP, rs2952768, is located closer to 
the METTL21A (FAM119A) gene than to the CREBl gene on 
chromosome 2 (Figures la and b). However, no significant 
association was found between rs2952768 and METTL21A mRNA 
expression, precluding the attribution of phenotypic alterations 
related to this SNP to the expression levels of the [\AETTL21A gene. 
Although the precise functions of METTL21A are poorly under- 
stood to date, a representative METTL, METTL11A, reportedly 
exhibited catalytic activity as a histone methyltransferase"^^ and 
chronic morphine treatment exhibited the acetylation and 
trimethylation of histones."^^ Although future studies are 
imperative, the action of opioids might be partially modulated 
by histone methylation via METTL functions. 
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With regard to our initial GWAS conducted as a consecutive 
three-stage analysis, the lowest combined P-value for the entire 
sample was P = 8.044 x 10"^ (Supplementary Table S2), which 
would have been deemed genome-wide nonsignificant if only a 
single-stage analysis was used to calculate conventional Bonfer- 
roni- or false discovery rate-corrected P-values for the total 
samples to determine statistical significance. However, 'significant' 
results obtained as conventionally corrected P-values will not 
always represent true associations, meaning that the results may 
not be necessarily replicated in other studies, and vice versa. For 
example, data from the National Human Genome Research 
Institute GWAS catalog (as of 31 January 2009), show 1321 entries 
of discovered associations with a P-value of < 10"^ but only 550 of 
these entries have a P-value of <5x10"^^^ which is a 
conventionally corrected conservative threshold for declaring a 
significant association in a GWAS.^^'^^ In both cases, truly potent 
candidate SNPs may be included in the outcome of the studies. 
Furthermore, GWASs in pharmacogenomics, such as this study, 
would tend not to yield 'significant' results obtained as 
conventionally corrected P-values compared with complex- 
disease GWASs^^ for several reasons. Among at least 16 different 
GWASs on drug response since the first was published in late 2007, 
less than half have shown genome-wide significance, although 
some potentially interesting associations that come close to 
significance have been detected in several of the studies in this 
category.^ ^ Altogether, these reports suggest that conventionally 
corrected P-values for the combined samples are not the only 
criteria to find true associations between SNPs and the phenotypes 
examined. The SNP we found, rs2952768, appears to be a 
promising SNP that is associated not only with opioid analgesic 
sensitivity in two independent surgical operations but also with 
several dependence-related traits in other subjects, prompting us 
to consider this SNP as the best candidate SNP known to date that 
is truly associated with human opioid sensitivity. 

Our compelling results suggest the possibility that the 
association observed in this study can be robustly generalized 
to various clinical and nonclinical scenarios, although this study is 
rather exploratory, and independent confirmation of the findings 
will be required in subsequent studies before various forms of 
practical clinical utilization of the prediction of opioid sensitivity 
based on this SNP can be applied. In conclusion, although the 
underlying mechanisms remain to be fully elucidated in future 
studies, our findings provide a novel step toward understanding 
individual differences in opioid sensitivity and stimulating future 
studies that can open new avenues for the personalized treatment 
of pain and drug dependence. 
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